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ABSTRACT 

Suggestive evidence has accumulated that intermediate mass black holes 
(IMBHs) exist in some globular clusters. As stars diffuse in the cluster, some 
will inevitable wander sufficiently close to the hole that they suffer tidal disrup- 
tion. An attractive feature of the IMBH hypothesis is its potential to disrupt 
not only solar-type stars but also compact white dwarf stars. Attention is given 
to the fate of white dwarfs that approach the hole close enough to be disrupted 
and compressed to such extent that explosive nuclear burning may be triggered. 
Precise modeling of the dynamics of the encounter coupled with a nuclear net- 
work allow for a realistic determination of the explosive energy release, and it is 
argued that ignition is a natural outcome for white dwarfs of all varieties passing 
well within the tidal radius. Although event rates are estimated to be signifi- 
cantly less than the rate of normal Type la supernovae, such encounters may be 
frequent enough in globular clusters harboring an IMBH to warrant a search for 
this new class of supernova. 

Subject headings: black hole physics - supernova: general - white dwarfs - glob- 
ular clusters: general 

1. Introduction 

The conjecture that there may be an intermediate mass black hole (IMBH) in some 
globular clusters has gained credence over the past few years. There is some dynamical 
evidence for a mass concentration within the central regions of some globular clusters. In 
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particular, the dynamics of stars in the inner regions of nearby clusters such as M15 and Gl 
suggest the presence of black holes with masses of the order of 10 3 and 10 4 M , respectively 
(IGerssen et al.ll2002l . 120031 ; iGebhardt et al.ll2002l . 120051 ) . This evidence remains rather contro- 
versial, partly because the velocity dispersion profiles can be reproduced without invoking the 
presence of an IMBH. Recently, somewhat more ambiguous evidence has arisen for the pres- 
ence of IMBHs in young star clusters, where ultraluminous, compact X-ray sources (ULXs) 
have been preferentially found to occur (IZezas et al.ll2002l ; iPooley k Rappaportll2006l ). Their 
high luminosities suggest that they are IMBHs rather than binaries co ntaining a normal stel 



lar mass black hole. Stud ies based on detailed iV-body simulations (IPortegies Zwart et al. 



2004 ; iPatruno et al.ll2006l ) of the evolution of a young star cluster in M82, the position of 
which coincides with an ULX, give further credence to this idea. While none of these argu- 
ments in itself carries enough weight to be convincing, the fact that they are so different in 
character does suggest that the existence of IMBHs in some young star clusters should be 
taken seriously. 

One would like some independent corroboration of the IMBH hypothesis, or, conversely, 
some way of ruling it out. Stellar disruption potentially offers such a test. As stellar orbits 
diffuse in phase space, it therefore seems inevitable that some will wander sufficiently close 
to the hole that they suffer tidal disruption. A star interacting with a massive black hole 
cannot be treated a s a point mass if it g ets so close to the hole that it becomes vulnerable to 
tidal deformations (IFrank k Reeslll976l ). The tidal radius, rx, defined as the distance within 
which a stars gets tidally disrupted, obviously has a value that depends on the type of star 
being considered. For a solar-type star it is 
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while for a degenerate white dwarf is roughly 
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Here r w d,g and M w d,o.6 denote the white dwarf radius and mass in units of 10 9 cm and O.6M , 
respectively. Note that for white dwarfs, the tidal radius is indeed inside the gravitational 
radius, r g ~ 1.5 x 10 M^ q cm, for black holes masses exceeding 
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For this reaso n, tidal disruption of whi t e dwarfs offers a unique d i agnostic for the presenc e 
of an IMBH flLuminet k PichorJ Il989l : Iwilson k Mathewsl booi baumgardt et all I2004T ). 
This paper explores the observational manifestations of such phenomena, with particular 
reference to white dwarfs that approach the hole close enough to be disrupted and severely 
compressed. 
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Numerical Methods 



Th e observat i onal consequences of stel lar disruption depend on what happens to the 
debris ( jReeall998l ; lEvans fc Kochanek!ll989l ). To quantify this, we have performed detailed 
three-dimensional, hydrodynamical calculations. The gas dynamics is coupled with a nuclear 
network to explore the effects of explosive nucleosynthesis during the strong compression 
phase. 

Hydrodynamics. The smoothed particle hydrodynamics method (SPH) is used here 
to solve the equations of hydrodynamics. Due to its Lagrangian nature SPH is perfectly 
suited to follow tidal disruption processes during which the corresponding geometry, densi- 
ties and time scales are chan ging y i olentl y. The SPH-formulation used in this study follows 
closely the one described in iBenzl (119901). Th e SPH equations derived from a Lagrangian 



( ISpringel Hernquist 112002 



Monaghanl 120021) yield di fferent symmetries in the particle in- 



dices, but a recent comparison (IRosswog fc Pried 120071 ) between these two sets of equations 
found only very mi nor differences. The forces from self-gravity are evaluated via a binary tree 
( iBenz et al.lll990bl ). the gravitation al forces from the central black hole are calculated using a 
Paczyhski-Wiita pseudo potential (iPaczvriski fc Wiitalll980T) . Details of how the singularity 
is treated numerically can be found in IRosswog I (120051 ). We have taken particular care to 
avoid artifacts from the use of artificial viscosity (AV). The so-called Balsara-switch (IBalsara 
19951 ) is implement ed here to avoid spurious s hear forces. In addition, time-dependent vis- 
cosity parameters (IMorris &: Monaghanl 119971 ) are used so that AV is only present where 



really needed. Details of the AV implementation and tests can be found in IRosswog et al. 
(feoooh . 



Equation of state. The HELMHOLTZ equation of state (EOS), developed by the Center 
for Astrophysical Thermonuclear Flashes at the University of Chicago is used. This EOS 
accepts an externally calculated chemical composition, facilitating the coupling to nuclear 
reaction networks. The electron-/positron equation of state makes no assumptions about the 
degree of degeneracy or relativity and the exact expressions are integrated numerically and 
tabul ated. The interpolation i n this table enforces thermodynamic consistency by construc- 
tion (jTimmes fc Swesty 1 12000| ) . The nuclei in the gas are treated as a Maxwell-Boltzmann 
gas, the photons as blackbody radiation. 

Nuclear burning. To account for the feedback onto the fl uid from nuclear transmutations 
we use a minimal nuclear reaction network developed by iHix et al.l (I1998I ). It couples a 
conventional a-network stretching from He to Si with a quasi-equilibrium-reduced a-network 
for temperatures in excess of 3 • 10 9 K. Although a set of only seven nuclear species is used, 
this network reproduces all burning stages from He- burning to nucle ar statistical equilibrium 
(NSE) accurately. For details and tests we refer to IHix et al.l (119981 ) . 
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3. Tidal Crushing and Ignition of White Dwarfs 

Several snapshots taken from our numerical simulations of a 0.2 M white dwarf ap- 
proaching a 1O 3 M black hole on a parabolic orbit with pericenter distance well within the 
tidal radius (r min = are shown in Figured! While falling inwards towards the hole, 

the star develops a quadrupole distortion which attains an amplitude of order unity by the 
time of disruption at r ~ rx- The resultant gravitational torque spins it up to a good fraction 
of its corotation angular velocity by the time it gets disrupted. The large surface velocities 
and the order unity tidal bulge combine to overcome the star's self gravity and lead to the 
disruption of the star. 

The behavior of a white dw arf passing well within the tidal radius exhibits special 



features (jLuminet fc Pichonlll989l ). As illustrated in Figured], the degenerate star is not only 
elongated along the orbital direction but is even more severely compressed into a prolate 
shape (i.e., a pancake aligned in the orbital plane). Each section of the star is squeezed 
through a point of maximum compression at a fixed point on the star's orbit. During this 
very short lived phase, the pressure grows sharply and the matter is raised to a higher 
adiabat. The distortions Ar w( j/r w d ~ 1 impart supersonic bulk flow velocities during the 
drastic compression of the stellar material. As a result, the temperatures increase sharply and 
trigger explosive burning (of He for the case shown in Figured]). The corresponding energy 
release can be very large (more than the star's self gravity) and it can build up sufficient 
pressure to significantly modify the dynamical evolution. Figure [2] shows the evolution of 
the compressed, and tidally disrupted white dwarf with and without the effects of nuclear 
burning. The pressure cannot build up sufficiently until after the white dwarf has endured a 
substantial volume contraction. This contraction is halted further by nuclear energy release. 
Just after the white dwarf attains its maximum compression, the density drops dramatically 
as a result of the subsequent expansion. 

During the brief but strong compression phase, the temperature reaches values beyond 
3 x 10 9 K, approaching but not quite reaching NSE. As a result, nuclei up to and beyond 
silicon are synthesized from the initially pure 4 He 0.2 M & white dwarf. The initial compo- 
sition was pure 4 He and the final mass fraction in iron-group nuclei is roughly 15%, located 
near the center (Figure [3]). This result should be taken as a modest underestimate, due to 
limitations in the seven species nuclear reaction network. Post-processing calculations, using 
a 300 isotope nuclear network over thermodynamic histories resulting from these calcula- 
tions, show significant nuclear flow above silicon, for helium- rich portions of the gas with 
peak temperatures above 2 x 10 9 K. As a result, heavier elements (like calcium, titanium and 
chromium) would be made, accompanied by a modest increase in the energy generation. It 
is thus safe to conclude that the white dwarf is tidally ignited and that a sizable mass of 
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iron-group nuclei is injected into the flow. Although the thermodynamical evolution and 
the nuclear energy release are sensitive to the initial stellar composition, here assumed to be 
4 He, we found that ignition is a natural outcome for white dwarfs of all masses (e.g. massive 
white dwarfs deprived of helium resulting from the evolution of a medium mass star) passing 
well within the tidal radius. For example, a C/O 0.6 M Q (1.2 M ) white dwarf approaching 
a 500 M black hole on a parabolic orbit with pericenter distance r min = r T /5 (r min = r T /2.6) 
ignites and, as a result, at least 0.17 (0.66) M Q of iron-group nuclei are synthesized in the 
flow. Thus, in the most favorable cases, the nuclear energy release may be comparable to 
that of typical type la supernovae. 

Although the explosion will increase the fraction of ejected debris, enough remains to 
be accreted on to the hole (Figure HJ). The returning gas does not immediately produce 
a flare of activity from the black hole. First material must enter quasi circular orbits and 
form an accretion torus jReesl l 19981 ; lEvans &: Kochaneklll989l ). Only then will viscous effects 



release enough binding energy to power a flare. The bound orbits are very eccentric, and 
the range of orbital periods is large. For white dwarfs, the orbital semi-major axis of the 

^ /o 2/3 

most tightly bound debris is a ~ 300M h3 r W d,9M wd06 r g , and the period is only t a ~ 
150(a/300r g ) 3 / 2 M h l^ 2 s. The simulation shows that the first material returns at a time 
< t a , with an infall rate roughly given by 1O 2 M yr _1 . Such high infall rates are expected 
to persist, relative steadily, for at least a few orbital periods, before all the highly bound 
material rains down. The vicinity of the hole would thereafter be fed solely by injection of 
the infalling matter at a rate that drops off roughly as t~ 5 ^ 3 for t > ~ 600 s. Once the 
torus is formed, it will evolve under the influence of viscosity, radiative cooling winds and 
time dependent mass inflow. A luminosity ~ -^Edd — 10 41 Mh 5 3 erg s _1 can therefore only 
be maintained for at most a year; thereafter the flare would rapidly fade. It is clear that 
most of the debris would be fed to the hole far more rapidly than it could be accepted if the 
radiative efficiency were high; much of the bound debris must either escape in a radiatively- 
driven outflow or be swallowed inefficiently The rise and the peak bolometric luminosity can 
be predicted with some confidence. However, the effective surface temperature (and thus 
the fraction of luminosity that emerges predominantly in the soft X-ray band) is harder to 
predict, as it depends on the size of the effective photosphere that shrouds the hole. 



4. Discussion 



Exactly how often a star gets close enough to the ce ntral hole to be tida lly disrupted 
depends on the stellar velocity distribution in star clusters (IFrank &: Reeslll976l ). In a simple 
case when the velocities are isotropic, the frequency with which a star enters the zone of 



- 6- 



vulnerability is 
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where is the star density in the star cluster nucleus and M^ 3 is the hole's mass in units 
of 10 3 M Q . This fiducial capture rate may need to be modified for a number of reasons. The 
actual rate could be lower than that given by equation (Tj0), even if the initial distribution 
were isotropic, if th e loss cone orbits were depleted faster than they could be replenished 
( IFrank fc Reeslll976l ). If the relaxation timescale were short enough to ensure replenishment 
of the loss cone, then the same diffusion processes could build up a cusp in the stellar 
distribution near the black hole. Capture or disruption of these stars could then augment the 
rate given by equation (J3]). None o f these complications is, however, likely to seriously modify 
the simple estimate given above (IBaumgardt et al.l I2004T ) . The rate given by equation (J4]), 
although highly simplified, agrees well with the fiducial r ates derived from detaile d N— body 
simulations of multimass star clusters containing IMBHs ( IBaumgardt et al.ll2004l ). However, 
in these simulations Baumgardt and collaborators found that IMBHs in star clusters disrupt 
mainly main sequence stars and giants, with white dwarfs accounting for only a few tenths 
of all disruptions. The probability of a pericentric distance within r m ; n is proportional to the 
first power of r min (see equation [1]). While complete disruption requires r min < rx, ignition, 
as calculated here, appears to be a natural outcome for white dwarfs passing well within the 
tidal radius: r min < r T /£ with £ > 3. The actual rate for white dwarf ignition would then 
be at least an order of magnitude lower than that given by equation (TJJ. 

It may seem, however, that even the modest rate of stellar disruptions given in equation 
flU) could have conspicuous consequences. The mass fraction that is ejected rather than swa l- 
lowed, though less spectacular than typical la supernovae (IHillebrandt &: Niemeyerl |2000| ). 
could nevertheless have distinctive observational signatures. The ejected material would 
be concentrated in a cone or fan close to the orbital plane (Figure [3]). The characteris- 
tic velocities in excess of ~ 2 x 10 4 km/s (Figure [3]), the kinetic output when a white 
dwarf is disrupted being ~ 4 x 10 50 erg s _1 . The light curve powered by radioactive decay, 
would indicate the synthesis of at least 0.1M w( j of 56 Ni before the spray of gas becomes 
translucent. With a rat e of disruption of white dwarfs £ ~ 1CT 8 yr -1 per globular cluster 
( IBaumgardt et al.l 12004 ). we expect the mean time interval between successive captures per 
galaxy to be ~ 10 5 (£/10~ 8 yr~ 1 ) _1 (/bh^V g c/10 3 ) _1 yr, where fhhN gc is the average number 
of globular clusters harboring an IMBH per galaxy and the number of globular clusters per 
host, iVgc, is found, for example, to vary between 5 x 10 2 and 7 x 10 3 from dwarf ellipticals 



to giant ellipticals (Ivan den Berghlll998l ; iBrodie fc Straderll2006l ). Event rates are then esti- 
mated to be significantly less than the rate of type la supernovae, which have an observed 



time interval between one type la explosion and the next of about 10 2 yr (ILeibundgutl 12001 
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Sharon et al.l 120071 ). but could be frequent enough to warrant a search for this new class of 
optical transient. 

The transient sky at faint magnitudes is poorly known, but there are major efforts un- 
der way that would in crease the discovery rate for type la supernovae from a few thousands 
(lAldering et al.l 120071 ) to about hundreds of thousands (IRiess fc Livio 1120061 ) per year. The 
transient itself should have several distinguishing characteristics. First, it should be seen in 
association with a globular cluster. Second, the explosion itself should be different, since the 
disrupted, degenerate stars should be, on average, lighter than those exploding as type la su- 
pernovae. Third, the spectra should exhibit large Doppler shifts, as the ejected debris would 
be expelled with speeds > 10 4 km/s and the optical light curve should be rather unique as a 
result of the radiating material being highly squeezed into the orbital plane (one thus expects 
different timescales for conversion of nuclear energy to observable luminosity when compared 
with normal type la events). Finally, the peculiar, underluminous thermonuclear explosion 
should be accompanied by a thermal transient (predominantly of soft X-ray emission) signal 
with peak intensity L ~ L-^dd = 10 41 M/ lj3 erg/s, fading within a year. Such transient signals, 
if detected, would be a compelling testimony for the presence of IMBHs in the centers of 
globular clusters. 
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Fig. 1. — A 0.2 M Q white dwarf (modeled with more than 4 x 10 6 SPH particles) approaching 
a 10 3 M Q black hole on a parabolic orbit with pericenter distance r m j n = is distorted, 

spun up during infall and then tidally disrupted. Temperatures (in units of 10 6 K) of 
the white dwarf-black hole encounter are shown. The panels in the upper row show cuts 
through the orbital (xy-) plane before and after passage through pericenter, as the white 
dwarf attains its maximum degree of compression. The panels in the lower row show the 
temperature distribution in the xz-pane (averaged along the y-direction) . 
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Fig. 2. — Evolution of the thermodynamical and mechanical properties of the central portion 
of the disrupted white dwarf. In this test calculation twenty central SPH-particles are chosen 
and their average temperature and volume (V = Y2i m i/Pi with rrii particle mass and pi par- 
ticle density) evolution is followed during passage through pericenter. As the star penetrates 
deeply within the tidal radius, the pressure cannot build up sufficiently until after the star 
has undergone an important volume contraction. The strong compression of the volume goes 
along with a sharp temperature increase. When burning is taken into account (red and green 
curves), the pressure terms grows very sharply, and, as a result, the temperature (averaged 
over the chosen set of particles) increases beyond 2.1 x 10 9 K. 
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Fig. 3. — The nature of the unbound debris at 110.6s. The panels show cuts through the 
orbital (xy-) plane well after passage through pericenter. The panels in the left columns 
show the mass fraction distributions of the synthesized group elements. The panels in the 
right column show the velocity distribution, (3 = v/c, and the specific energy, e, (in code 
units). 
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Fig. 4. — Influence of the nuclear burning on the ejected matter: differential mass distribu- 
tions in specific energy for the 0.2 M & white dwarf debris. The resultant explosive energy 
severely decreases the amount of material bound to the hole from 50% to 35% of the initial 
mass of the white dwarf. 



